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Abstract. Holes in a Mott insulator are represented by spinless fermions in the fermion-boson 
^ , model introduced by Edwards. Although the physically interesting regime is for low to moderate 

$_^ ' fermion density the model has interesting properties over the whole density range. It has 

jy^ ' previously been studied at half-filling in the one-dimensional (ID) case by numerical methods, 

in particular exact diagonalization and density matrix renormalization group (DMRG). In 
C^ ' the present study the one-particle Green's function is calculated analytically by means of a 

decoupling scheme for the equations of motion, valid for arbitrary density in ID, 2D and 3D 
with fairly large boson energy and zero boson relaxation parameter. The Green's function is used 
' O ' to compute some ground state properties, and the one-fermion spectral function, for fermion 

C , densities n=0.1, 0.5 and 0.9 in the ID case. The results are generally in good agreement 

with numerical results obtained by DMRG and dynamical DMRG and new light is shed on the 
nature of the ground state at different fillings. The Green's function approximation is sufficiently 
successful in ID to justify future application to the 2D and 3D cases. 
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1. Introduction 

Hubbard-like models provide a paradigm for a large class of strongly correlated systems. 



t — 

f— ^ ■ general form for the Hubbard Hamiltonian is 



^Hu = T + U = - ^ tpabipfeCra + U^ 11^+12^- (1) 

r,p,a,b r 



and the system is strongly correlated when the repulsive on-site interaction U is considerably 

5^ I larger than the hopping parameters t. Here Cra destroys an electron in state a on lattice site r 

and cl^ , creates an electron in state 6 on a nearest-neighbour site r -|- p. The state indices a, b 

are summed over two states denoted by -|- and — , and the occupation numbers nr± = c^^Cri. In 
general we consider bipartite lattices in one, two and three dimensions. The standard Hubbard 
model [IJ corresponds to the case of a single orbital on each site with states it corresponding to 
spin ±1/2 and with hopping parameter 



'^pab 



to^ab ■ (2) 



This model with a 2D square lattice is frequently used to describe the copper-oxygen plane of 
high-Tc systems such as doped La2Cu04, where the orbital corresponds to a Cu dx2_y2 orbital. 



A related model on a 2D square lattice in the xy plane describes a ferromagnetic system in which 
the on-site states all have the same spin and the states + and — correspond to d^2_^2 and d22_y2 
t2g orbitals respectively. This "t2g model" describes ferromagnetic planes in Sr2V04 and also in 
fluorides such as K2CUF4 and Cs2AgF4, where in this case a crystal field converts eg orbitals into 
an effective t2g system [2H6] . In the t2g model one considers only the dominant hopping processes 
in which hopping between + states occurs along the x axis and hopping between — states occurs 
along the y axis. Thus in the t2g model 

tpab = tQ5ah\p • GqI , (3) 

where e+ and e_ are unit vectors along the x and y directions respectively and p is a unit vector 
in the nearest-neighbour direction p. 

In the large-?/ limit we can approximate the Hamiltonian ?^hu by an effective Hamiltonian 
T^eff which is defined to act only within a subspace where there is no double occupation of any 
site. T-LeS takes the form 

nes = T-T^/U, (4) 

where 7i is that part of T which has matrix elements linking states of no double occupation 
with states having a single doubly-occupied site. If three-site terms are neglected it is straight- 
forward to write the second term of Eq. (j3|) in terms of spin (pseudospin for the t2g model) 
operators S^ = (nr+ — rir-)/2, S^ = S^ ± Sr = cl^Crz^ and number operator Ur = nr+ + Ur-. 
Then, for the Hubbard model, TieS becomes 

/it— J = /t-J +2*^ / ^ i'-'r'-'r+p + 2 V'-'r ^r+p + '-'r ^r+p) ~ 4^r^r+pJ 
r,P 

= Tt-J + i J ^(Sr • Sr+p - InrUr+p) , (5) 



r,P 



where J = 4tQ/C/ and 7t-j contains the hopping parameters of Eq. ([2]) [7j. For the t2g model 
the corresponding Hamiltonian is 

^t2g = 7t2g + 4 J / ^('S'r'S'r+p — J^-r^^r+p) , (6) 

r,P 

where Ttjg contains the hopping parameters of Eq. ^ [2j. 7^t-j is known as the t-J model and if 
the transverse exchange terms in Eq. 1^ are omitted we have the t-J^ model. This differs from 
"Htjg where -|- and — spins only hop along the x and y axes respectively. For the case of one 
electron per site, the undoped case for the oxides being modelled, T drops out and the system is 
a Mott insulator. In the Hubbard case the insulator is a Heisenberg antiferromagnet and in the 
t2g case it is an Ising antiferromagnet with alternating orbital order in the ground state. The 
transverse exchange terms in Eq. ([5]) are not present in Eq. ([6|) owing to the directional hopping 
of the orbitals in the t2g case. 

When a hole is introduced into the lattice its motion disturbs the spin or orbital order of 
the ground state. At each hop through the ordered lattice the hole leaves a spin, or orbital, 
deviation at the site it vacates. Thus as the hole hops through the lattice the energy of the 
system increases linearly and the hole is bound to its starting-point. This is known as the string 
effect and it exists in dimensions higher than one [8]. In ID the spin deviation created by a 
hole at its first hop, from some initial position, increases the energy but subsequent hops in the 
same direction merely shift the ordered spin configuration by a lattice spacing without further 
increase in energy. There is therefore no string effect in the ID t-J or t-J^ models. Furthermore 
in the 2D or 3D t-J models the string effect is relaxed by the S^S~j^p terms in Eq. ([5]) which 



exchange spin directions between lattice sites and can lead to a healing of the spin deviations 
created by the hole. Thus the hole can move at a speed determined by the healing rate, which 
leads to a quasiparticle band of width proportional to the exchange parameter J [9] . This way 
of relaxing the string effect is not possible in the t-J^ or the ^t2g model owing to the absence of 
the S^ S^j^p terms. In these models a relaxation mechanism can be introduced by including in 
the Hamiltonian a term of the form 

Y,{St + S-) = 2Y,S^, (7) 

r r 

so that the spin part of Eq. ([6]) becomes a transverse-field Ising antiferromagnet. In the physical 
case of the t2g model this transverse field corresponds to an on-site crystal field which mixes 
the two t2g orbitals just as the transverse magnetic field mixes + and — spins. In both the 
t-J^ and the %i^ model inclusion of three-site terms mentioned earlier also relaxes the string 
effect [2] but we shall not include them here. A hole in the 2D t-J^ model can also propagate by 
means of a Trugman path |10] which consists of 6 hops around a 4-site square plaquette. The 
hole moves to a next-nearest neighbour site and leaves the antiferromagnetic spin arrangement 
undisturbed. This cannot occur in the 'Ht2g, model owing to the directional hopping [2]. We 
shall here concentrate on the t-J model and the t-J^ model in a transverse field although there 
should be no difficulty in treating the directional hopping of the ^t2g niodel. 

Our reference state has one electron on each site with spins (or pseudospins) ordered as in an 
antiferromagnet. We follow Martinez and Horsch [IJL] in introducing a spinless fermion operator 
/r which creates a hole in the reference state at site r and a boson operator 6j which creates a 
spin reversal on site r. Thus 6j = S~ for site r on the + spin sublattice and hi = S^ for r on the 
— sublattice. Clearly each nearest-neighbour hop of the hole through the ordered lattice reverses 
the spin on the site vacated by the hole. Thus the hopping operator 7t-j may be written as 

7I-J = -toE(/r+p/r^J + H-C.)- (8) 

rp 

For the t2g model the corresponding hopping operator Ttjg is 

7^2^ = -io Y.^fUpfA + H.c.) \p ■ ej . (9) 

rp 

Since S^ = ±(2 — fcrfcr) for r on the it spin sublattice we may write the Ising part of the exchange 
term, in the presence of a hole, as 

- i J ^^(1 - hr){l - 26t6r)(l - 25|:+/r+p)(l - /ir+p) = ^^bl^il - h^+p - bJ+Zr+p) , (10) 
rp rp 

where h^ = fifr and an irrelevant constant term has been dropped. Here we have used hrblb,. = 
since the hole and a spin deviation cannot occupy the same site. A spin deviation created on 
site r when the hole vacates that site will in general have one neighbouring site r + p occupied 
by the hole and another by a spin deviation created when the hole arrived on site r. Thus the 
expression in Eq. (jlOp may be written as 



ij(z-2)^6t6r, (11) 

r 

where z is the number of nearest neighbours, z = 2 in ID so that the above term is zero, which 
is consistent with the absence of the string effect in the ID t-J^ model. Apart from an additive 



constant the t-J Hamiltonian Eq. ([5]) becomes 

^t-J = -to Y.(fr+pfrbl + H.C.) + ljY,(blbU, + 6r&r+p) + ^(^ " 2) J^ bl^ ■ (12) 

rp rp r 

Note that the derivation of this Hamiltonian in Ref. [llj gives additional constraints which we 
will not discuss here. 

The relaxation of the string effect occurs when bosons, created by fermion hopping, are 
spontaneously destroyed in nearest neighbour pairs by the terms &r^r+p- Edwards [12] introduced 
a simplified model Hamiltonian of the form 

-^Ed = -to ^ifl+pfA + H.C.) - X^ibl + 6r) +^oYl blK + NX^ojo (13) 

rp r r 

in which boson relaxation terms 6r^r+p are replaced by the simpler linear ones br- A unitary 
transformation "^li^Ed = e^Ede"*^, S = (X/ojo) X^r(^r — ^r), resulting in b,. ^ br + A/cjq, yields 

-flEd = -^ Yl /r+p/r - to E(/r+p/r4 + H.C.) + 6^0 J] blb^ . (14) 

rp rp r 

Thus the second term in "HEd is eliminated in favour of the first term in T^Ed which introduces a 
coherent hopping channel in addition to the original incoherent one. In the ground state |0) of 
the Hamiltonian ?^Ed) in the absence of fermions, there are no bosons, so that 6r|0) = 0. Thus 
for the ground state |0) of ?^Ed we have {br — A/li;o)|0) = and hence (0|6r5r|0) = A^/wq. 

Prom the above discussion it is clear that in 2D the Edwards model corresponds, within a 
certain range of parameters, to an underlying t-J^ type of model with Hamiltonian 

-toYl 4+pCr. + ^^0 Yl SrS^p - 2A j; S^r , (15) 

rpcr rp r 

since 25^ = Sr + Sr = bl + br. This corresponds to a doped antiferromagnetic Ising model in a 
transverse magnetic field, but this correspondence is only valid for low hole density (low fermion 
density in the ^Ed model) and A <C Wo- In this case, since KOjS'rlO)! = ^ — {0\brbr\0) = ^ — X"^ /loq, 
the holes move in a background medium which is close to a saturated antiferromagnet. For 
X > ojo the transverse-field Ising model no longer exhibits antiferromagnetic order [13]. In ID 
the correspondence breaks down completely since if A = and wq 7^ the 'HEd model exhibits the 
string effect which is absent in the ID t-J^ model. In fact the ID T-LEd niodel has an interesting 
similarity to the 2D t-J^ model. Even for A = the string effect is relaxed by an analogue 
of the Trugman path [T^. This is again a 6-step process in which the fermion propagates to 
a next-nearest neighbour site with the background medium left undisturbed, i.e. no bosons 
excited. During the process the fermion excites three bosons which are subsequently destroyed. 
Unfortunately the model T^Ed does not represent precisely any physical system that we know 
of. The underlying Ising exchange is characteristic of the t2g model for vanadates and fluorides 
but the directional hopping of that model is not included, although this could be remedied in 
the 2D model. However if the model is considered over the whole X/to, too/ to space it is found, 
even in the ID case, to exhibit a surprising number of different physical regimes reminiscent 
of some found in realistic strongly-correlated electron systems and strongly-coupled electron- 
phonon systems. This was first demonstrated for the case of a single fermion (Nj = 1) at 
temperature T=0 by the method of variational Lanczos diagonalisation [TJ]. A great advantage 
of the model is that the simple treatment of the background medium, in terms of local bosons. 



makes it possible to obtain essentially exact results in the thermodynamic limit, at least in the 
ID case. More recently the half-filled case, with Nj = N/2 where N is the number of lattice 
sites, has been investigated |15H17j . The ground-state phase diagram has been mapped out in 
the whole X — loq plane, using a density matrix renormalisation group (DMRG) technique |18pi9j . 
A quantum phase transition between a metallic Tomonaga-Luttinger liquid and an insulating 
charge density wave (CDW) was shown to exist. 

It is desirable to complement these numerical results with some more analytical approaches. 
Nearly all analytical work on t-J-like or polaronic models is confined to the case of a single 
fermion. However very recently the projective renormalisation method [2U] has been applied to 
the half-filled ID 7^ Ed niodel in a study of CDW formation at T=0 [21]. Our initial analytical 
work on the ^Ed niodel was confined to the one-fermion case in ID at T=0 with the additional 
restriction A = |22j . The latter restriction means that the fermion can only propagate 
as a coherent quasiparticle by the Trugman-like process discussed above. This requires the 
coexistence of at least three bosons during the particle's motion and in Ref. [22] the one-fermion 
Green's function was calculated within a 3-boson approximation. The spectral functions agree 
well with exact numerical results for wq/^o ^ 1- In this paper we extend the analytical approach 
for A = to finite fermion density and to higher dimension, 2D and 3D. This is achieved within 
a 2-boson approximation and comparison with exact numerical results in ID shows that the 
range of validity is thereby reduced to wq/^o ^ 3. No Trugman-like processes exist in the 2- 
boson approximation but they are unimportant to understand CDW formation at finite density. 
Relaxation of the string effect does not depend on them because bosons created by one fermion 
can be destroyed by other fermions. Near half-filling a CDW state is found, in good agreement 
with ID numerical results, and the CDW transition temperature can be calculated formally. 
However, within the 2-boson approximation, the transition is mean-field-like, with no short- 
range order in the disordered state. This is clearly not realistic in ID but finite-temperature 
calculations should certainly be relevant in 3D. In this paper we concentrate on the calculation 
of spectral functions in ID at T=0 with various fermion densities and compare with some new 
DMRG results. 

In Sec. [2] we determine the one-fermion Green's function for A = within the 2-boson 
approximation in the cases of one particle (Nj = 1) and one hole (Nj = N — 1). The method 
used is different from that previously employed for Nf = 1 [22] and provides a wave-function as 
a byproduct. In Sec. [3]we study the hierarchy of equations of motion for the Green's function at 
arbitrary density and devise a decoupling which leads to the correct Green's function in the 2- 
boson approximation for Nf = 1 and Nf = N — 1. This 2-boson Green's function is valid for any 
bipartite lattice in ID, 2D or 3D and allows for two distinct sublattices so that the CDW state 
can be investigated. In Sec. H] spectral functions and ground-state properties calculated from 
the 2-boson Green's function are compared with those calculated numerically by the DMRG 
method. In Sec. Ejwe draw conclusions and consider the outlook for future work. 

2. Limiting cases of the Green's function ^vithin the 2-boson approximation 

We now derive expressions for the one-fermion Green's function in the special cases of a single 
fermion {Nf = 1) and a single hole {Nf = N — 1). For these single-particle cases the Green's 
function may be calculated by direct solution of the Schrodinger equation, which also yields the 
wave- function. 

2.1. The case Nf = 1 

For A = the Hamiltonian Eq. p3p . written in k-space, takes the form 



^kk 



where 



/k 



^Ee'^-Vr 






7(k) = -toJ]e^- 



kp 



(17) 



The p-summation is over z nearest neighbours in a ID, 2D or 3D bipartite lattice. In the 2-boson 
approximation the wave-function for a single fermion is of the form 



^i 



/.^ + E«(qi)/k-q.^l + EE«(qi'q2)/, 



^ 6^ At 

k-qi-q2^qi"q2 



vac) 



(18) 



qi 



qi q2 



where |vac) is the vacuum state and a(qi,q2) = a(q2,qi). On substituting this in the 
Schrodinger equation ^^k = -E^k) and multiplying on the left by (vac|, (vac|6q//k-q' and 
(vac|6q/ 6q/ /k-q' -q' we obtain equations of the form 



qi 



1 2 ^-^ 

=7(k-qi) + a(qi)(wo--E) + ^= > a(qi, q2)7(k - qi - q2 

N VN q^ 



7(k-qi -q2)[a(qi) + a(q2)] + 2a(qi, q2)(2tJo - E) = . 



On solving Eq. (pT|) for a(qi,q2), and substituting in Eq. ([20]) . we find 



o(qi) 



loq — E 



zti 



2iOo - E 



' 7(k-qi) + ^E-''"'''"'''^^'"^'''^ 



N 



q2 



2ujo- E 



(19) 

(20) 
(21) 

(22) 



It is easily shown that for a bipartite lattice the last term in this equation vanishes when 
a(q) oc 7(k — q). Hence a solution for a(qi) is obtained by omitting the last term and, on 
substituting this solution in Eq. (fT9]) . we find 



E- 



zti 



ujQ — E 



zt\ 



0. 



(23) 



2a;o-_E 



The solutions of this equation are the energies of single- fermion eigenstates, which are the poles of 
the one- fermion Green's function Gk(-E'). The left-hand side of Eq. ()23p is G^ (E). The absence 
of k-dependence in this expression shows that the single-fermion eigenstates are localized. This 
is due to the string effect which is not relaxed in the 2-boson approximation. 



2.2. The case Nj = N - 1 

The motion of a single hole in the present model is quite different from that of a single particle. 
To avoid confusion it should be stressed that the hole discussed here does not correspond to a 
hole in a t-J-like model, the latter being represented by a fermion in the present model. Clearly 
when the hole considered here hops to a neighbouring site a boson is created on the arrival site, 
not on the departure site as in the motion of a single particle. This boson can be destroyed 
immediately when the hole makes a further hop. There is therefore no string effect and the hole 
propagates easily. To find the Green's function in this case we use the Schrodinger equation as 
in the case of a single particle. The wave-function for the hole is of the form 



$1 



/k + y^c(qi)/, 



qi 



k+qi 



^qi+EE^(^l''l2)/k+qi 



b^ b^ 
+q2%'^q2 



\F), 



(24) 



qi q2 



where \F) is the state with every site occupied by a fermion and c(qi,q2) = c(q2,qi). The 
equations corresponding to Eqs. (fT9]) -([2T ]) in the previous case are 

-E--^7(k)j;c(qi)=0, (25) 



-^7(k) + c(qi)(^o -E)- -^7(k + qi) Yl c(qi, q2) = , (26) 



'c(qi)7(k + qi) + c(q2)7(k + qs)] + 2c(qi, q2)(2wo - E) = . (27) 



From Eqs. ()26p -()27 p we find that for a bipartite lattice 

. X 1 7(k) , , 

cqi =^ ^Vt^ -y. 28 

^.o-E-hM^ 

Hence, from Eq. ([25]). 

-E-[7(k)]2ly ^ r . M2 =0- (29) 

The solutions of this equation are the energies of single-hole eigenstates and, as expected, they 
depend on the wave- vector k, being functions of [7(k)]'^. This type of k dependence arises from 
the fact that the hole propagates through the lattice, leaving behind no excited bosons, by 
means of double hops, as discussed at the beginning of this section. As in the one-particle case 
we deduce that for one hole {Nf = N — 1) the one-fermion Green's function is given by 

[0.(E)l-=£-Mk)plE— ^^H^- PO) 

q Wo + ^ - 2^J^;+E 



The sign of E has been changed since in Eq. ()29p the energy refers to a hole state. 

3. The Green's function at finite fermion density 

In this section we study the hierarchy of equations of motion of the one-fermion Green's function 
and find a decoupling which is consistent with the results derived in Sec. [2] for the limiting cases 
of low (Nf = 1) and high (Nf = N — 1) fermion density. We allow for two distinct sublattices 
with different occupation so that the CDW state can be investigated. 

The Fourier transform of the one-fermion retarded Green's function is defined by [23] 



CUE) = iih-jl)) = -i / rfie(t)([/k(i),/k]+)e'^*, (31) 



where /k(i) = e'^*/ke '^* and 9{t) is the unit step function. We may write /k, defined by 
Eq. pT|) . as a sum of two sublattice components. Thus 

/k = ^(/ki + /k2), (32) 

where 

^' = \/IE«*-''^. (33) 



Here the suinination is over sites r/ which belong to sublattice / {I = 1,2). It follows from 
Eqs. dSID and ([32]) that 

G'k(i^) = ^EE^k'"(^)' (34) 

1=1 m=l 

where 

G^ = ((/kz; /L)) = I E «''^-^''"''"^Gr,s,„ , (35) 

rjSm 

and 

Gr,s„ = ((/rp/L)). (36) 

The equation of motion for the Green's function {{A; B)) is |23j 

E{{A-B)) = {[A,B]+) + {{[A,n\-B)). (37) 

Hence, noting that % is given by Eq. ()13p with A = 0, we find 

p p 

where 

^ri+p,s,n ^ iifri+pbri', fs^)) , (39) 

4^;,t: = ((/r3U;/L)>. (40) 

Here r^ is a general site on the sublattice I which complements the sublattice I (I = 2, 1 for 
/ = 1, 2 respectively), r^ +p is a particular site on the sublattice 7 but it is necessary to define / 
more generally in order to close the equations of motion. However we first calculate the second 
term on the right of Eq. ()38p . 

The equation of motion for H is 

{E - CJo)i^r'+p,s„ = Sri+p,sm{bvt) - to E^^^+P-^n+p'/r, 5 /!„)) 

(41) 

-to ^{{fr,+p+p'iKi+p + bl+p+p')K^-jL)) ■ 
P' 

The last term of this equation is produced by a process in which a fermion hops from r/ + p to 
ri+ p + p' with either the creation of a boson on the vacated site r^ +p or destruction of a boson 
on the arrival site r; + p + p'. The latter boson must have been created by another fermion 
so that the term involving 6j , „ , „/&r; corresponds to a dynamical interaction between fermions. 
This should only be included if we consider two-fermion interactions consistently which goes 
beyond the effective Hartree-Fock treatment introduced below. We therefore neglect this term. 
In the second term on the right of Eq. ()4ip we retain only the part involving fermions on two 
sites, thus taking p' = p, and then make a Hartree-Fock type of approximation. Thus the second 
term becomes 

- to(((l - nr,+p)/r,; /D) ^ -toil - (nr,+p))Gr,,s„ , (42) 

where nY^ = fJifr,- Also we write the average fermion occupation for sites on sublattice / as 
ni = (nri) so that (nr,+p) = nj. Following this discussion Eq. ([^T]l becomes 

{E - a;o)^;;+p,s„ = Sr,+p,s^ (br,) -toil- nj)Gr„s^ - to E "^n W,s™ ' (43) 



where 

'^;;:;;>,s. = ((/r,+p+p'^.+p^. ; /L» • (44) 

In the equation of motion for J the terms arising from [fn+p+p' ^T~L\ involve three boson operators 
and we neglect them. The terms arising from [fcrj+p^rp^] niay be treated in the Hartree-Fock- 
like way used to obtain Eq. (fl2]) . Hence 



{E - 2u.o)j;;+;;';.,s„ = 5r,+p+p',.^K+phr,) - hi - nz)i/;;+p,s„ . (45) 

We approximate (&r;+p^r;) by 5^6^, where 6/ = {hn)- Furthermore it is shown below that 6/ = 
for the present model, with A = 0. Hence {hn+phn) may be neglected. On substituting for J in 
Eq. (j43p . using Eq. (j45p . and summing over p, we find 

E Hl\+p,s^ = C{E, m) [bl Y, Sr,+p,s^ - zt,{l - nj)Gr„s^] , (46) 

p p 

where 

The second term on the right of Eq. (j38p has therefore been determined and to find the third 
term we must consider the equation of motion of the Green's function I defined by Eq. (j40p . 
This takes the form 

[E + U^)n!.+P^ = 5rj,s^{bl^+p)+toY.{{fr,fUpfr,+p+p'-jL)) 

-to E(( Vp'('5r-,r,+p + 4,+p^r, + ^j^^+p' 4, +p ) ; /sl» • 
P' 

We now make similar approximations to those used to obtain Eq. (|43p . Thus we make an 
effective Hartree-Fock approximation to the second term on the right of Eq. ()48p . retaining it 
only when rj = r^ + p, and we neglect the b^ ^p^r- term in the third term. Hence 

P' P' 

where 

<S':- ' = ^^ Vp'6U,4,+p; /D) • (50) 

We treat the equation of motion of if in a similar way to that of J. A slight difference is 
that instead of neglecting all the terms arising from [fr-+p' , H] we retain one which leads to a 

factor br-+pib^^_j^ , . This equals [l + br__|_./&r-+p'] which we approximate by 1, neglecting the boson 

occupation number. Hence, neglecting a (6^6^) correlation function as before, we find 

{E + 2.,)K^2:' = -^on. E ^:^+P",s. • (51) 

p" 

Thus, from Eqs. ([M]) and dSU, 

{E + a;o)/-t: = '^r.s^: - tonA,,+p E ^r^+p'.. + ^f|- E E <;,,».. • (^2) 

p' p' p" 



To solve this equation for / we introduce the Fourier transform 



3 

^fe = (I) ' E ^riZ e'[('^+'')---k-'-'i-(r<+P)l (53) 

so that 

(^ + -0 - ^^^mS^) ll^ = /f V'^l™^* + /|-n(k)GL™ . (54) 

By taking the Fourier transform of Eq. ()38p , and using Eq. ()46p , we find 

EG'^ = di^ - C{E, ni)[ba{k)&i^ - ztlil - n^)G(r] + 7(1^)/^ , (55) 

where 

rySm q 

By combining the last three equations we obtain 

G^^[E-ztl{l-nj)CiE,ni) - [j{k)]\D{E,ni)] = 6im - [C(^,nz)6; + D(^, k,nz)6f](^^„7(k) , 

(57) 
where 

D(E, k, ni) = riTT^T , (58) 

E + u}o- -B+2ZJr 

and 

DiE,ni) = ^^DiE,q,ni). (59) 

q 

It remains to explain why, as indicated following Eq. ()45p . 6/ = 0. This may be shown using 

■ ' i-oo (60) 

= ^ E J_jHrUp,r,+f,{E + ir?) - i^rVr.+p(£^ " ^dE . 

Eqs. (j46p and (j57p may be used to show that the above expression is a linear combination of 
bi and 62. Thus 6/, b^ (/ = 1, 2) satisfy a set of linear homogeneous equations so that in general 
bi = 0. This result is consistent with the symmetry of the transverse-field t-J^ model (Eq. ()15p ) 
which underlies the present model for low fermion density. When A = 0, as assumed here, the 
transverse magnetic field vanishes so that by symmetry the expected value of the transverse spin 
moment (S^) = {bl + bj.)/2 = 0. The properties bi = and the characteristic periodicity in k 
are shown in the Appendix to be generally true for the model with A = for all values of wq/^o- 
As soon as A 7^ we shall in general have 6^ 7^ so that, from Eq. ()57p . the k dependence of the 
Green's function will involve 7(k) and not only [7(k)]^ as is the case for A = 0. 
The final result for the Green's function in the case A = is 

G'^{E) = 5i^[E - ztlil - nj)C{E,ni) - ^ikfnjDiE,ni)]-' (61) 

where C and D are given by Eqs. ()47p and ()59p respectively. Also, from Eq. ()34p . 

GUE) = ^[G^i{E) + G,,2{E)] , (62) 



where Gkj = Gj^- In the homogeneous case, where there is no CDW, ni = 722 = n, where 
n = Nf/N is the fermion density, and Gki = Gk2 = G^. For n = we recover the A'^^ = 1 result 
for G^ , given by the left-hand side of Eq. (j23p , and for n = 1 we recover the Nf = N — 1 result 
given by Eq. ()30p . Thus Gk is correct in these two limits within the 2-boson approximation. 
The factor Sim in Eq. ()6ip shows that within the present approximation fermions propagate, by 
means of double hops, within a single sublattice. The dimension of the bipartite lattice (1, 2 
or 3) enters only through the band energy 7(k) [Eq. ()17p ] and the function D, which may be 
written as 

D{E, m) = [ ^f^].^^^ / 7 , (63) 

where NQ{'y) is the density of states per site for the band energy 7(k). 

To complete the present formulation we require the equations which determine the chemical 
potential fi and the self-consistent sublattice densities ni,n2 for a given fermion density n. The 
density of states per site on sublattice I is given by 

Ni{E) = -^Yl I™ ^k« {E + ir?) (64) 



iVvr 

k 



and 



711= f Ni{E)f{E,fi)dE, (65) 

where f{E,fi) = [e^'^"'^) -|- 1]~^ is the Fermi function with chemical potential /x and /3 = 
{kBT)~^. The chemical potential and the CDW order parameter P are determined by 

^{ni+n2)=n, ^(ni - 77-2) = nP . (66) 

If the system is ordered (P 7^ 0) at T = the CDW temperature Tc, where P — )• 0, can be 
calculated. In the present mean- field- like approximation this will not be the true Tc but a higher 
temperature where short-range order substantially disappears. The calculated Tc should be a 
reasonable approximation to the true value in 3D. In general, if quantities such as the spectral 
functions discussed below are calculated assuming P = 0, the results will be valid in the high 
temperature limit where there is truly no short-range order. 

The spectral function for states projected onto the sublattice I is given by 

Siik,E) = --lmGM{E + iv) (67) 

vr 

and the total spectral function is 

Sik,E) = -^luiGk{E + iri) = hsiik,E) + S2{k,E)]. (68) 

vr z 

These spectral functions may be used to calculate the occupation number n(k) = (/j[/k) and 
the related quantity d(k) = (/k/k+Q/2)' where Q is a basis vector of the reciprocal lattice, /k 
is given by Eq. ([32]) and it follows from Eq. ([33]) that /k+Q/2 = (/ki — fk2)/V^- Hence we find 

n(k) = 1 r [5i(k, E) + S2(k, E)]dE (69) 

^ J-00 

and 

d{k) = \ r [5i(k, E) - S2(k, E)]dE . (70) 



These quantities satisfy the sum rules 



n 



^^d{k)=nP. 



(71) 



In the next section we report numerical results for all the above quantities, based on the 
approximate Green's function of Eq. (|6ip . in the ID case. Many of the results are compared 
with those of the DMRG method, both to assess the validity of the present approximation and 
sometimes to throw new light on the DMRG results. The application of the DMRG method to 
the present model has been described previously |16yi7j. 

4. Numerical results for ground-state and spectral properties in ID 

In ID the band energy 7(k) = — 2cosA;, with the lattice constant taken as 1 and the unit of 
energy taken as to- The quantity D which appears in the Green's function (Eq. (j61|) ) can then 
be evaluated analytically using Eqs. (j58p and (j59|) . The q-summation is most conveniently 
performed as a contour integral around the unit circle. The result is 



DiE,ni) = DirOiaf ^ bf) ~ iDiiOib 
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The density of states, given by Eq. ()64p . can be evaluated similarly with the result 

Ni{E) = Na{E) + Ni2iE) , 
where 
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Ni^{E) = -e{a 
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2njDii{E) 



(72) 
(73) 

(74) 

(75) 

(76) 

(77) 



(78) 
(79) 

(80) 



The results presented below are for the case ojq = 3. By comparison with the DMRG results 
for the ground state properties n(fc) and d{k), and with the dynamical DMRG ^24j (DDMRG) 
results for S{k,E), it is found that this boson energy is large enough for many results of the 
present analytic approximation to be quite accurate. 



4.1- High and low fermion density 

The existence of a two-sub-lattice CDW state in the present model is well-established at half- 
filling (n = 0.5) for sufficiently large wo/^o |15ftl7t l21j. However the general results of the 
Appendix imply that a two-sublattice CDW state is always a possibility, whatever the density 
n. We therefore used Eqs. ()65p . ()66p and ()75p to search for such states even in the high and low 
density cases of n = 0.9 and 0.1. In the case n = 0.9 with wo/to = 3 we find a self-consistent 
CDW state with order parameter P = 0.063. The fiow of the iterative procedure to determine P 
indicates that this, not the uniform density P = state, is the stable ground state. In the DMRG 
calculations convergence to the CDW solution is improved by use of suitable external fields at the 
boundaries with open boundary conditions. The spectral function S{k, E) calculated from the 
Green's function using Eq. (j68p is compared with the DDMRG results in Fig.[TJ The quasiparticle 
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Figure 1. Zero-temperature single-particle spectral function S{k, E) at n = 0.9 for loq = 3 from 
Eq. (|68p (left panel) compared with the DDMRG result (right panel) obtained for a finite system 
with A^ = 20 sites using open boundary conditions (OBC). The inset shows the photoemission 
spectra (PES) and inverse photoemission spectra (IPES) near the Fermi point. E is measured 
with respect to the Fermi energy Ep (all energies are given in units of to)- 



peaks for the Green's function decoupling scheme results (left panels) are delta-functions in the 
limit r/ — )• 0, but to make them visible we have taken rj = 0.05; the same value has been taken in 
the DDMRG data (right panel). The following main features are in good agreement: the general 
shape and width of the quasiparticle band crossing the Fermi level, the dispersive peaks just 
below E — Ep = —4 which vanish for k = 7r/2, the weak flat band at E — Ep = 2 in the left panel 
and 3 in the right one. The splitting of the quasiparticle peaks due to CDW order is clearly 
visible in the left panel. This splitting is not clearly resolved in the main right panel but the inset 
for k = 7r/21 shows a splitting between a peak below Ep in the photoemission spectrum (PES) 
and one above Ep in the inverse photoemission spectrum (IPES). The absence of the splitting 
for states further from Ep in the DDMRG data is presumably due to finite lifetime broadening 
processes which are not included in the Green's function approximation. It should be noticed 
that the quasiparticle states at the Fermi level are from only one of the split subbands near k = Q 
and vr, the other subband being fully occupied. This means that the Fermi wavevectors are at 
vr/lO, TT — 7r/10 so that the hole pockets near k = Q and tt have the correct Fermi surface 'volume', 
this being twice what it would have been in the absence of the CDW {P = 0). Discontinuities 
at these wavevectors are clearly seen in Fig. 3 where n{k) and d{k), calculated from Eqs. ()69p 
and ()70p . are plotted. The agreement with the DMRG data, also plotted, is excellent. In the 
Green's function calculations the discontinuities are not perfectly sharp owing to the use of a 
finite r/ (=0.0005) near Ep. There are two contributions to n{k) and d{k), one arising from the 
coherent quasiparticle bands and the other from spectral weight further below the Fermi level. In 
the present case d{k) arises almost entirely from the coherent contribution and the quasiparticle 



states have slightly more weight on the minority sublattice, hence the negative value of d{k) 
over most of the zone. Near k = and vr, however, only the subband associated with the 
majority sublattice is occupied, hence the strong positive contribution. Correlations giving rise 
to Luttinger-liquid behaviour in ID are beyond the present Green's function approach. It should 
be noted that at n = 0.9 the CDW state is metallic, whereas at n = 0.5 it is an insulator, as 
discussed in Sec. 14.21 

In the case n = 0.1 with ujo/to = 3 we find a self-consistent CDW state with order parameter 
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Figure 2. Single-particle spectral function S{k,E) at n = 0.1 for A = and ojq = 3 from 
Eq. dMD. 
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Figure 3. Bloch state occupation numbers n{k) (left) and related correlation function d{k) 
(right) for n = 0.9 (upper panels) and n = 0.1 (lower panels) at A = 0, wq = 3. Analytical 
results (lines) are compared to numerical DMRG data (symbols). 



P = 0.61. The discrepancy in order of magnitude between this value and the much smaher one 
for n = 0.9 is due to the consistent use of fermion density in the definition of P (Eq. ()66p ). rather 
than changing to hole density for the n = 0.9 case. For n = 0.1 it has proved difficult to converge 
to a CDW solution in DMRG. This may indicate that the Green's function approximation is 
failing in this low density case. If so, the nature of the ground state is unclear. Nevertheless 
in Fig. [2] we show results for the spectral function and in Fig. [3] we plot n{k) and d{k). The 
bottom of the narrow quasiparticle band near Ep is at k = 7r/2 and only the majority subband 
is occupied. Hence d{k) is strongly positive near k = tt/2 and almost equal to n{k). Over most 
of the zone the only contribution is from spectral weight further below the Fermi level; this 
vanishes at k = 7r/2 and makes a negative contribution to d{k) owing to more weight residing on 
the minority sublattice. The discontinuities in n{k) and d{k) at the expected values 7r/2 it vr/lO 
are clearly seen. 



4.2. The CDW state at half- filling 

We now consider a self-consistent CDW state for the half- filled band with n = 0.5 and A = 0, 
ojQ = 3. As before the Fermi energy Ep (or chemical potential /i at T > 0) and order parameter 
P are determined by Eqs. (f66|) and (f65]) . with the sublattice densities of states given by Eq. ([75]) . 
For the above parameters we find /i = — 0.211, P = 0.765. Ep lies in a gap as is appropriate for 
an insulator. 

The spectral functions calculated from the Green's function method and the DDMRG are 
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Figure 4. Single-particle spectral function S{k, E) at n = 0.5 for X = 0, ojq = 3 from Eq. ([68 
(left panel) compared with the DDMRG result (right panel). 
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Figure 5. Density of states D{E) near the Fermi level for A = at half filling from Eq. ([75 
(left panel) and corresponding quasiparticle band dispersion E{k) (right panel). 
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Figure 6. Bloch state occupation number n{k) (left panel) and related correlation function 
d{k) (right panel) for the half-filled band case with X = 0, ojq = 3. 

compared in Fig. HI The agreement is generally good. The Fermi level lies in a gap between 
two quasiparticle bands, a broad upper unoccupied one and a narrower lower occupied one. The 
main discrepancy is that in the DDMRG case the lower quasiparticle band is extremely narrow 
whereas in the Green's function method it has a significant width. This shortcoming of the 
Green's function approximation is discussed thoroughly in Sec. 14.31 Fig. [5] shows the density of 
states near the Fermi level, projected onto the majority and minority sublattices, calculated from 
Eq. ()64p . It is remarkable that states in the occupied quasiparticle band are entirely confined 
to the majority sublattice, whereas those in the unoccupied quasiparticle band reside entirely 
on the minority sublattice. The dispersion curves E{k) of the quasiparticle bands, obtained 
by plotting the loci of quasiparticle peaks in the spectral function, are shown in FiglSj The 
approximate Green's function therefore predicts an indirect gap, with the top of the occupied 
quasiparticle band at A; = 0, vr and the bottom of the unoccupied band at A; = 7r/2. 
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Figure 7. Single-particle spectral function (left panel) and quasiparticle band (right panel) in 
the high-temperature limit P = for the half-filled band case with A = 0, wq = 3. 



In Fig. [6] left and right panels we plot curves for n{k) and d{k), respectively, calculated from 
Eqs. ([69]) and ([70]) . together with points obtained by the DMRG method for finite systems of 
different sizes. The agreement is very striking, particularly since in DMRG these quantities 
are calculated directly from the ground state whereas in the Green's function method they 
are obtained as integrals over the spectral function. It should be noted that n(7r/2) = d{7r/2) 
because for k = 7r/2 all the spectral weight of the occupied states resides in the quasiparticle 
peak which is entirely based on the majority sublattice. From Eqs. ([69]) and ([70]l it is clear that 
the weight in this peak is 2n(7r/2) ~ 0.88. 



Finally in Fig. [7] we plot the spectral function (left panel) and quasiparticle energy (right 
panel) in the high-temperature limit, where all short-range order has disappeared (P = 0). The 
Fermi wave- vectors are close to 7r/4 and 37r/4 with the correct Fermi surface 'volume'. This is 
the state from which the CDW evolves as the temperature is lowered. Clearly the situation is 
quite different from the usual one in which a CDW evolves because of nesting between Fermi 
wave-vectors at k = ib7r/2. The origin of this difference is that for A = only next-nearest 
neighbour hopping, resulting in vr-periodicity in fc-space, occurs even in the disordered state. 

4-3. Discussion 

The most notable difference between the Green's function and DDMRG results lies in the width 
of the occupied quasiparticle band for n = 0.5. A very narrow band is also found in earlier 
calculations for coq = 2 using exact diagonalisation where there is also no tendency for the top 
of the band to be at fc = 0,7r rather than tt/2 [15]. The band seems to be almost as narrow 
as one would have from the Trugman-like six-step process in a perfect CDW {P = 1). This is 
quite surprising since for ujq = 2 the densities on the two sublattices are about 0.8 and 0.2 [16] 
which is far from a perfect CDW. Satellites below the main quasiparticle peak are suggestive of 
finite-size effects in the systems considered with 12 and 16 sites. 

To gain more insight into the Green's function approximation we may consider the limit of 
large ujq where all the weight is concentrated in the two quasiparticle bands. From Eq. ()6ip it 
follows that these bands are given by 

,,, ±2P + (1 ± P) cos(2k) , , 

E{k) = ^^ '- ^^, (81) 

LJQ 

where the upper signs (+) correspond to the upper band and the lower signs (— ) to the lower 
band. The widths 2(1 it P)/ujq of these bands can be understood by inspection of the following 
processes, where a bullet represents a vacant site, a star represents a boson and a circle represents 
a fermion: 

I • 00)^I0 *0)^I0 0-) 
IOO-)^IO *0)^l • OO) 

In the upper line the central site is on the majority sublattice. The diagram shows how a fermion 
added to the upper band, on the minority sublattice, can hop by a two-step process across an 
occupied majority site. The probability of the majority site being occupied is (1 -|- P)/2 so an 
estimate of the width of the upper band is 2(1 -|- P)/ujq, as in Eq. (|81|) . In the lower line the 
central site is on the minority sublattice. Clearly a hole created on the majority sublattice can 
hop by a two-step process across an occupied minority site. The probability of the minority 
site being occupied is (1 — P)/2 so an estimate of the width of the lower quasiparticle band is 
2(1— P)/a;o as in Eq. (|8ip . The probability argument used here is equivalent to the Hartree-Fock- 
like approximation made in deriving the Green's function. This approximation is clearly failing 
in the present situation, since the width of the lower quasiparticle band is much larger than that 
given by the DDMRG. To expose the cause of this failure we used DMRG to calculate the three- 
site correlation functions (nm-in-m^-m+i) for a 10-site system with periodic boundary conditions 
and A = 0, ujo/tQ = 3 as usual. The correlation function takes two values 6.46062 x 10~^ and 
3.90826 X 10~^, depending on whether m is an even or odd site corresponding to the majority and 
minority sublattice, respectively. The simple probability argument would give respective values 
P(l - P)2 = 0.04225, (1 - P)p2 = 0.1375 when P = 0.765. If a site on the minority sublattice is 
occupied it is clear that the probability of both neighbouring sites being occupied is only about 
0.039. Thus one of these majority sublattice sites is very likely to be occupied by a hole. Hence 
there is a strong tendency for the minority sublattice fermions and the majority sublattice holes, 
equal in number, to form bound pairs. We conclude that for n = 0.5 important correlations 



exist in the ground state which are not included in the Hartree-Fock-hke approximation of the 
Green's function decoupUng. 

5. Conclusions and outlook 

Doped Mott insulators remain at the forefront of physics research largely because of their 
relevance to high temperature superconductivity. In the cuprate superconductors holes in the 
copper oxide planes move in a background of antiferromagnetic order. In other systems of 
interest the background is one of alternating orbital order. The progress of a hole through such 
a background is hindered by the string effect. This effect has been known for a long time [8] 
but the string picture is central to much recent work (e.g. |25fl28j ). The Edwards fermion-boson 
model considered in this paper was introduced to describe this effect in the simplest possible 
way [12]. The spinless fermions correspond physically to the holes in the Mott insulator. The 
ordered background does not appear explicitly in the model; the essence of the string effect 
actually relies only on the existence of substantial short-range order without the necessity of 
true long-range order. Clearly the physical interest lies mainly in 2D [29j, but so far most 
calculations for the Edwards model have been made in ID |14ftl71[2T| l22). In this paper we 
describe an analytical approximation to the one-fermion Green's function which is valid in ID, 
2D and 3D. Its main limitations are that the boson energy should be fairly large (wq/^o > 2) and 
that string relaxation is neglected (A = 0). The principal objective of this paper is to test the 
accuracy of the Green's function method, within its expected domain of validity, by comparing 
with numerical results obtained in ID by the DMRG and DDRMG methods. 

This comparison has been made in detail for wq = 3 and for various fermion densities. For 
the half-filled band case (n = 0.5) excellent agreement is obtained for ground state properties. 
These include the CDW order parameter, the Bloch state occupation number n{k) and a related 
quantity d{k) associated with the CDW. There is also generally good agreement for the one- 
fermion spectral function although the Green's function method predicts much too wide an 
occupied quasiparticle band. In Sec. 14.31 the reason for this discrepancy is traced to missing 
correlations in the Green's function approximation. The Green's function method predicts that 
the CDW state evolves from a high-temperature disordered state with Fermi wave-vectors close 
to 7r/4 and 37r/4, which is quite different from the usual case where a CDW evolves because of 
nesting between Fermi wave- vectors at A; = ib7r/2. The origin of the difference is that only next- 
nearest neighbour hopping, resulting in vr-periodicity in A;-space, occurs even in the disordered 
state. 

There is excellent agreement between the two methods for the dilute hole (n = 0.9) case with 
the usual parameters A = 0, wo/^o = 3. Somewhat surprisingly, there is also a two-sublattice 
CDW state which in this case is metallic. In the Appendix it is shown how the appearance of 
two-sublattice CDW states is related to a symmetry property of the model with A = 0. The 
comparison between the two methods proved to be more difficult in the dilute fermion (n = 0.1) 
case. The Green's function method again predicts a metallic CDW state but the DMRG fails 
to confirm this. The true nature of the ground state in this case remains unclear. 

It may be concluded that the rather simple Green's function approximation derived here 
is sufficiently successful in ID, which is probably the least favourable case, to envisage future 
applications to the 2D t — J^ model and, with a slight extension of the model, to the t2g model 
of alternating orbital order. The situation of physical interest will be low to moderate fermion 
density, the fermions corresponding to holes in the relevant Mott insulator. 
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Appendix 

Through the particular form of the fermion-boson hopping term at A = (cf. Eq. ^ 
the transfer of a fermion beween neighbouring lattice sites coincides with a change of the 
number of bosons by one. As a consequence we may think of the bosons as tracking the 
motion of the fermions. This picture becomes exact through the idenfication of a conserved 
quantity. Let us define operators N^ and N? which count the number of fermions on the 
A- or B-sites of a bipartite lattice, and similar operators N^ and N^ for bosons. Then 
Nfh = Nf — N? + 2{N^ — N^) commutes with the Hamiltonian HEd, when A = 0, and is 
therefore a conserved quantity. The eigenvalues of Nfh can be used to classify the eigenstates of 
the Hamiltonian. We note that Nfb is not conserved for A 7^ 0. 

The existence of Nfb has two major consequences. First, fermion operators such as 
fnAfniB change Nfb by 2, so that in any eigenstate of the Hamiltonian the expectation value 
ifnAfmB) ~ ^ ^^^ arbitrary sites nA and mB on the respective sublattices. Similarly (6„) = 0. 
Hence quantities such as the spectral function ^(k, E) or the momentum distribution n(k) have 
the periodicity in k of the reciprocal lattice of the real-space A or B sublattice, e.g. vr-periodicity 
in the ID case. Second, it implies that those eigenvalues of the Hamiltonian corresponding to 
Nfb 7^ are degenerate. This follows because the translation operator T commutes with H, 
but changes the sign of -/V/b- The energy eigenvalue of an eigenstate \^/J) with {^|J\Nfb\1l^) 7^ 
must therefore be (at least) two- fold degenerate, since the state T\ip) belongs to the same energy 
but differs from \^p) due to the change of Nfb- This degeneracy corresponds to a breaking of 
translational symmetry, as for a two-sublattice CDW state. Clearly if the ground state at A = 
is not such a CDW state it must have Nfb = 0. 
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